IDC  tii[  nof’ 

MEASUREMENT  OF  WALL  SHEAR  AND  WALL  PRESSURE 

DOWNSTREAM  OF  A 

HONEYCOMB  BOUNDARY  LAYER  MANIPULATOR 

by 

JAMES  CHRISTIAN  MOLLER  and  PATRICK  LEEHEY 


Report  No.  97457-3 
April,  1989 

This  research  was  carried  out  under  the  Underwater  Acoustics  Program 
Code  1 1250A  of  the  Office  of  Naval  Research  under  contract  number  i 

N0(X)14-86-K-0183 .  j 

Approved  for  public  release  ;  distribution  unlimited  13 

1 

Acoustics  and  Vibration  Laboratory 


O) 

CM 

00 

O 

CM 

< 

I 

D 

< 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  D.i.  Enitrtai 


REPORT  DOCUMENTATION  PAGE 


I.  report  number  ;  GOVT  accession  no 

Acoustics  &  Vibration  Laboratory 
^ - 1  No.  97457-3  _  _ _ 


I 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  catalog  number 


4.  TITLE  (fyd  Subtitle) 

5  TYPE  or  REPORT  A  PERIOD  COVERED 

Measurement  of  Wall  Shear  and  Wall  Pressure 

Downstream  of  a  Honeycomb  Boundary  Layer 

Technical 

6  performing  org.  report  number 

Manipulator 

7.  AUTMORr»> 

6  CONTRACT  OR  GRANT  NUMBER^^J 

James  Christian  Holler  and  'Patrick  Leehey 

N00014-86-R-0183 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 
Code  1125  OA  Ballston  Tower  #1 

800  N.  Quincy  St.,  Ai.Tington.  VA  22215-5000 _ 

M.  monitoring  agency  NAMC  a  AOORESSf/f  dlttetcr,!  C  on  f  ^  o  / 1 . 


10.  PROGRAM  ELEMENT.  PROjeCT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


ONR  Code  11250A 


12.  REPORT  DATE 

April  1989 _ 

13  NUMBER  or  PAGES 

_ 73 _ 

'5  security  CV.ASS.  (o!  rmport) 


Unci assif ied 


<$/»  oecl  Assin  cation  / down  grading 
SCh  EOULE 


‘6.  OJSTRIBUTION  statement  (0l  thia  Hapoft) 


Approved  for  public  release;  distribution  unlimited. 


»7.  O1STRI0UTION  statement  (ot  th0  4bHrmel  ainltfad  In  Dtock  20.  II  dllltrent  from  Raport) 


KEY  WORDS  (Conlinu*  on  rmvorao  sido  it  noc  *  *  Mmry  and  Identity  by  biock  numbar) 

Wall  shear  stress 

Wall  pressure  fluctuations 

Turbulent  boundary  layer  manipulation 


20.  ABSTRACT  (Comlnuo  on  rovergo  gido  i(  nncfusory  o/n/  uli'iMily  by  nuni6o»^ 


See  page  2  of  this  report, 


_ _ l^NCI.ASSIFIED 

SCCuniTY  CLASSIFICATI  ON  OF  THIS  P  A  f.  F  CR>..n 

89  5  31  038 


DD  \  JAN  73  1473  EDITION  OF  I  NOV  6'  IS  OBSOLETE 


MEASUREMENT  OF  WALL  SHEAR  AND  WALL  PRESSURE 

DOWNSTREAM  OF  A 

HONEYCOMB  BOUNDARY  LAYER  MANIPULATOR 


by 


JAMES  CHRISTIAN  MOLLER  and  PATRICK  LEEHEY 
Report  No.  97457-3 
April,  1989 

This  research  was  carried  out  under  the  Underwater  Acoustics  Program 
Code  1 1250A  of  the  Office  of  Naval  Research  under  contract  number 

N00014-86-K-0183  . 


Approved  for  public  release  ;  distribution  unlimited 


ACOUSTICS  AND  VIBRATION  LABORATORY 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


MEASUREMENT  OF  WALL  SHEAR  AND  WALL  PRESSURE 
DOWNSTREAM  OF  A 

HONEYCOMB  BOUNDARY  LAYER  MANIPULATOR 


by 

James  Christian  Moller  and  Patrick  Lechey 


ABSTRACT 


A  boundary  layer  manipulator  in  the  shape  of  a  honeycomb 
placed  in  a  turbulent  boundary  layer  has  been  shown  experimentally 
to  reduce  local  wall  shear  downstream  by  as  much  as  38%  and  wall 
pressure  by  36%.  The  zone  of  measureable  shear  and  pressure  reduc¬ 
tion  extends  for  at  least  300  pre-manipulator  displacement  thicknesses 
downstream.  Pressure  spectra  of  the  manipulated  flow  scales  better 
with  inner-outer  variables  than  does  the  unmanipulated  flow.  Con¬ 
vective  pressure  speeds  are  reduced  by  as  much  as  ^proximately 
30%  for  Strouhal  numbers above  0.3.  Below  = 

0.3  convective  phase  speeds  axe  not  noticeably  changed.  This  indi¬ 
cates  pressure  disturbances  for  <o8*/f/^  <  0.3  may  originate  from  the 
inner  portion  of  the  boundary  layer  since  the  manipulator  scales  are 
far  larger  than  the  viscous  sublayer.  Coherence  appears  to  be  reduced 
in  the  stream  wise  direction  and  increased  in  the  spanwise  direction. 

All  tests  were  done  in  the  same  facility  as  in  Blake'  {1}  using 
the  same  microphones  and  wall  insert.  Unlike  his  work  and  similar  to 
the  work  of  other  more  recent  investigations'*  {2,3,4}>.  peaks  were 
found  in  plots  of  coherence  vs.  phase  in  the  vicinity  of 
®5*/Ug,=0.3.  , 
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1,  INTRODUCTION 

1.1.  Background 

The  use  of  devices  to  alter  and  control  the  structure  of  turbulent  boundary  layers,  most  noteably 
the  Large  Eddy  Breakup  Device  (or  LEBU).  have  been  studied  by  several  researchers  in  recent  years. 
The  primary  goal  of  this  research  has  been  to  investigate  the  viability  of  the  LEBU  as  a  drag  reduction 
device.  While  an  understanding  of  the  flow  behavior  around  a  LEBU  is  reasonably  well  understood, 
results  have  been  mixed  and  seem  to  differ  because  of  differences  in  lest  facilities  [5).  While  some  have 
achieved  an  overall  (plate  plus  device)  drag  reduction,  others  with  seemingly  similar  experimental  rigs 
have  not  gotten  the  same  results.  The  idea  of  manipulating  the  boundary  layer  is  still  a  viable  one,  how¬ 
ever,  from  the  standpoint  of  wall  pressure  reduction  because  manipulators  do  create  local  (but  not 
overall)  reductions  in  wall  shear  for  many  boundary  layers  downstream.  Knowledge  of  boundary  layer 
wall  pressure  spectra  and  the  means  of  controlling  it  have  several  aero-  and  hydro-acoustic  applications. 
Perhaps  the  foremost  areas  are  for  the  prediction  and  cor.?,  ot  of  flow-induced  vibration  and  the  control  of 
boundary  layer  self  noise  over  a  hydrophone  array.  Beeler  [6]  has  reported  net  reductions  of  wall  pres¬ 
sure  of  12.5%  downstream  of  a  tandem  LEBU.  The  gains  to  be  had  from  boundary  layer  control  can  be 
significant.  For  example,  detection  probability  by  a  submarine  can  conceivably  be  increased  40%  by  a 
10  db  increase  in  signal-to-noise  ratio  [7]. 

The  linkage  between  wall  shear  and  wall  shear  pressure  is  as  follows.  It  is  hypothesized  wall 
pressure  scales  with  mean  wall  shear  stress.  If  we  begin  with  the  Navier-Stokes  and  continuity  equations, 

+  vV2(7 

Di  p 

and 

1.  Blake,  W.  K.,  (1970),  "Turbulent  boundary  layer  wall-presture  Ruclualions  on  smooth  and  rough  walls",  J.  Fluid 
Mechanics,  vol.  44,  part  4,  pp.  637.660. 

2.  Farabee,  Theodore  M.,  (1986),  "An  Experimental  Investigation  of  Wall  Pressure  Fluctuations  Beneath  Non-Equilibrium 
Flows",  D.  W.  Taylor  Naval  Ship  R  A  D  Center,  Kept.  No.  DTNSRDC-86/047,  pg.  126.  (May  1986) 

3.  Brooks,  T.F.  and  Hodgson,  T.M..  (1981),  "Trailing  Edge  Noise  Prediction  From  Measures  Surface  Pressures",  J.  Sound 
and  Vibration,  vol.  78,  No.  1,  pp.  69-1 17. 

4.  Narayan,  N.  and  Plunkeu,  R.,  (1985),  "Pressure  Cross  Spectra  in  turbulent  Boundary  Layers  in  Water",  ASME  Winter 
Annual  Mtg.,  .Miami  Beach,  Ha..  .Nov.  17-22,  NCA-Vol.  I.  pp.  97-103. 

5.  B^c^^'n,  Eugene  F.,  "European  Mtg.  on  Turbulent  Drag  Reduaion",  ESN  41  -2,  1987,  pg.  89. 
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vi7  =  o 

Taking  the  divergence  of  both  sides  gives 

V  ■  +  V  ■  (f7  •  V)f/  =  -±  V2n  +  vV  •  V2U 

^  p 

where 

V  •  (yiH)  =  V  ■  [V{V  ■  U)  -  Vx(Vxf7)] 

The  first  term  on  the  right  hand  side  is  zero  due  to  continuity  and  the  second  is  zero  by  the  rules  of  vec¬ 
tor  algebra.  Also  by  continuity  and  vector  algebra, 

=  =  0 

Which  leaves  the  Poisson  equation, 

V  (f7  V)£7-H  .lV2p  =0 
or 

After  rewriting  pressure  and  velocity  in  terms  of  mean  and  fluctuating  quantities  and  requiring  the  flow 
Geld  to  have  a  single  nonzero  mean  vector  component,  the  equation  to  describe  fluctuating  pressure  in  a 
boundary  layer  is 

V2p-  =  -p[2 -  Ui'uj)  ] 
oy  ux  ^Xf  3xj 

The  first  term  inside  the  bracket  is  known  as  the  mean  shear-turbulence  interaction  term  and  has  been 
experimentally  found  to  be  the  dominant  source  of  pressure  fluctuations.  In  it  is  a  link  to  wall  shear 
because  mean  wall  shear  is  linearly  related  to  the  mean  velocity  gradient  at  the  wall.  The  second  part  of 
the  mean  shear  turbulence  term,  d\/dx,  was  beyond  the  scope  of  this  set  of  experiments  but  would  need 
to  be  measured  in  order  to  get  a  complete  idea  of  the  scaling  of  the  wall  pressure. 

The  results  of  some  previous  researchers  also  indicate  that  the  wall  pressure  spectrum  scales  with 
mean  wall  shear  stress,  especially  for  frequencies  corresponding  to  the  smaller  structures  in  the  flow  [2]. 

6.  Beeler.  G.  B.,  (1986),  "Turbulent  Boundery-Liyer  Will  Pressure  Ructuations  Downstream  of  Tandem  LEBU",  AIAA 
Journal,  vol.  24,  no,  4,  pp.  689-691,  (April  1986). 

7.  Cox,  Albert  W..  Sonar  and  Underwater  Sound  .  Lexington  Books,  Lexington,  MA  (1974),  pg.  1 17. 
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1.2.  Objectives 

The  objectives  of  this  work  were  1)  to  find  the  manipulating  device  which  created  the  most  rapid 
and  largest  decrease  in  local  friction  coefficient  downstream  of  the  device,  2)  to  measure  and  describe 
the  changes  in  the  boundary  layer  across  and  downstream  of  the  device,  3)  to  make  point  measurements 
of  mean  wall  shear  and  wall  pressure  spectra,  4)  to  verify  that  pressure  spectra  scale  with  wall  shear  in 
the  non-equilibrium  flow,  5)  to  make  streamwise  and  spanwise  cross-spectral  measurements  in  the  non¬ 
equilibrium  flow  for  the  sake  of  describing  the  convective  behavior  and  structure  of  the  flow. 


2.  EXPERIMENT  DESCRIPTION 

2.1.  The  Device 

A  honeycomb  placed  inside  the  boundary  layer  was  chosen  as  the  device  best  suited  for  the  goal 
of  rapid  downstream  local  friction  reduction.  A  survey  of  the  local  firiction  coefficient  profiles  under 
manipulated  boundary  layers  reported  by  other  researchers  [8,9,10,11]  was  done  as  a  means  of  selecting 
likely  candidates  for  the  manipulator  desired  for  these  experiments.  Bushnell  makes  brief  mention  of  a 
honeycomb  placed  inside  the  turbulent  boundary  layer  offering  as  much  as  50%  reduction  in  wall  shear. 

Any  device  which  alters  the  large  scale  structure  of  the  boundary  layer  through  the  use  of  stream- 
wise  elements  in  the  outer  portion  of  the  boundary  layer  can  be  thought  of  as  having  two  effects;  "wake" 
effects  and  "plate"  effects.  The  wake  which  the  element  sheds  must  have  vorticity  in  both  senses  and 
that  which  is  opposite  to  the  sense  of  the  boundary  layer  forms  a  barrier  to  the  penetration  of  momentum 
from  structures  in  outer  portion  of  the  boundary  layer  to  the  inner.  The  presence  of  an  element  in  the 
boundary  layer  creates  a  new  no-slip  condition  and  thereby  serves  to  suppress  fluctuating  velocities.  The 
presence  of  these  two  conditions  ought  to  create  a  relatively  quiescent,  "re-laminarizcd"  region 

8.  Hefner.  J.  N.,  Weinstein,  L.  M,.  &  Bushnell,  D.  M.,  (1980),  "Ijirge-Eddy  Breakup  Scheme  for  Turhulcnt  Viscous  Drag 
Reduction",  Progress  in  Astronautics  and  Aeronautics ,  vol.  72,  pp.  110-127. 

9.  Anders,  J.  B.,  Hefner,  J.  N.,  and  Bushnell,  D,  M..  (1984).  "Performance  of  l.arge  Etldy  Breakup  Devices  at  Post- 
Transitional  Reynolds  Numbers",  AIAA  Paper  84-0345,  (Jan.  1984). 

10.  Corke,  T.  C.,  Guezennec,  Y.,  and  Nagib,  H.  M..  (1979).  "Modification  in  Drag  of  Turbulent  Boundary  l.ayers  Resulting 
from  Manipulation  of  Large-Scales  Structures",  Symposium  on  Viscous  Drag  Reduaion,  Dallas,  Tex.,  Nov.  7-8. 

1 1.  Erei,  D..  and  Thomann,  H.,  (1980),  "Direct  Measurement  of  Skin  Friction  in  a  Turbulent  boundary  I.aycr  with  a  Strong 
Adverse  Pressure  Gradient",  J.  Fluid  Mechanics,  vol.  101,  part  1.  pp.  79-95. 
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downstream  of  tliC  device.  A  honeycomb  may  be  less  effective  at  shedding  counter-vorticity  and  more 
effective  at  breaking  up  large  structures  than  the  LEBU  because  of  the  honeycomb’s  vertical  elements. 

A  honeycomb  essentially  introduces  three  new  length  scales  into  the  description  of  the  manipulated 
boundary  layer;  the  comb  height,  diameter,  and  streamwise  length.  Figures  la  and  lb  are  a  sketch  and  a 
photograph  of  the  device.  In  this  experiment,  all  three  were  kept  constant  and  the  boundary  layer  thick¬ 
ness  incident  upon  the  manipulator  was  varied.  A  fourth  variable  introduced  by  the  jM’esence  of  the 
manipulator  is  the  downstream  distance  from  it  to  the  point  of  measurement.  The  honeycomb  dimensions 
were  chosen  based  on  the  boundary  layer  thickness  and  the  corresponding  large  eddy  sizes.  The 
"natural"  boundary  layer  thickness  on  the  walls  of  the  duct  was  estimated  through  use  of  the  code  given 
by  Cebeci  [12],  The  comb  height  was  chosen  to  be  no  higher  than  the  boundary  layer  thickness  in  order 
to  not  introduce  any  irrotationality  into  the  essentially  inviscid  freestream.  The  hole  size  {d*=4O0)  was 
chosen  so  that  the  largest  eddies  would  be  most  affected  (i.e.  the  comb  height  to  hole  size  ratio  would  be 
less  than  ten).  This  is  much  larger  than  the  Taylor  microscale  for  the  Reynolds  numbers  of  the  experi¬ 
ment  After  Lumley  [13], 

The  boundary  layer  was  calculated  to  be  ^proximately  one  inch  in  thickness  before  entering  the 
device.  The  comb  used  was  made  from  aluminum,  had  a  mean  hole  size  of  1/4",  and  a  web  thickness  of 
.004".  The  honeycomb’s  streamwise  length  varied  with  distance  from  the  wall  with  the  narrowest  part 
being  in  the  outer  part  of  the  layer  and  the  widest  part  being  at  the  wall.  This  was  done  to  avoid  creating 
enough  blockage  to  separate  the  flow  and  create  a  quasi-recirculating  region  downstream  of  the  device. 
The  order  of  the  pressure  drop  in  the  outermost  parts  of  the  comb  can  be  estimated  by. 


4c.  = 


ndcIcVi^^ 


8v/. 


0.0008. 


This  profile  also  serves  to  be  screen-like  (tending  to  decrease  turbulent  length  scale)  in  the  outer  portions 
of  the  boundary  layer  where  eddies  are  relatively  large  and  velocity  fluctuations  relatively  low  and  tube- 


12.  Cebeci,  T.  and  Bradshaw,  P.,  Momentum  Transfer  In  Boundary  Layers  ,  Hemisphere  Pub.  Corp.,  Washington, 
(1977),  p.  194. 

13.  Tennekes,  H.,  and  Lumley,  J.  L,  A  First  Course  In  Turbulence  ,  M.l.T.  Press,  Cambridge,  Ma.  (1981),  pg.  67. 
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like  (lending  to  suppress  fluctuating  velocities)  in  the  inner  portions  of  the  layer  where  length  scales  are 
small  and  fluctuating  velocities  are  large. 

2.2.  Wind  Tunnel 

All  tests  were  carried  out  in  the  M.I.T.  Acoustic  and  Vibration  Laboratory  low-turbulence,  low- 
noise  wind  tunnel.  It  was  constructed  and  described  by  Hanson  [14]  and  is  shown  in  Figure  2.  The 
tunnel’s  blower  motor  has  recently  been  replaced  by  a  General  Electric  Model  DC-300  speed  controlled 
motor.  Drift  of  the  blower  speed  from  the  desired  speed  set  point  no  longer  occurs  with  the  new  motor. 
The  experiments  described  here  were  run  with  the  tunnel  in  the  closed  duct  configuration.  The  duct  is 
15'’xl5"  in  cross-sectional  area  and  approximately  13*  in  length.  The  maximum  core  flow  velocity 
through  the  duct  is  approximately  SO  meters  per  second. 

Four  tunnel  nominal  mean  flow  speeds  were  chosen  for  most  of  the  tests;  two  the  same  as 
Farabee’s  [2]  operating  points  (16  m/s  and  28  m/s),  one  between  the  two  (25  m/s),  and  the  fourth  at  a 
higher  speed  (33  m/s).  They  shall  hereafter  be  referred  to  as  speeds  1, 2,  3,  and  4,  respectively. 

The  honeycomb  manipulator  was  fastened  with  fine  wire  to  the  wall  of  the  tunnel  on  the  side 
which  was  upstream  of  the  shear-  and  pressure-measuring  insert  which  will  be  described  later.  The  data 
are  presented  for  various  streamwise  locations  downstream  of  the  honeycomb  and  this  was  accomplished 
by  moving  the  honeycomb  to  different  stations  along  the  tunnel  wall  rather  than  moving  the  insert  The 
honeycomb  was  placed  at  three  stations  upstream  of  the  insert  which  shall  be  referred  to  as  stations  1 ,  2, 
and  3  and  correspond  to  approximately  17.S,  SO,  and  80  pre-manipulator  displacement  thicknesses 
separation  between  honeycomb  and  measuring  station. 

All  cracks  between  tunnel  sections  and  at  the  instrumentation  insert  were  puttied  and  sanded 
smooth  with  600  grit  sandpaper.  The  tunnel  wall  downstream  of  the  comb  was  sanded  as  well. 


t4.  Hanson,  C.  E.,  (1969),  "The  Design  and  Constniction  of  a  Low  Noise,  Low  Turbulence  Wind  Tunnel",  Acoustics  and 
Vibration  Lab  M.I.T.  Rept.  no.  7961 1  -1 . 
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23.  Pressure  Measurements 

All  pressure  measurements  were  made  on  the  wall  of  the  duct  at  Reynolds  numbers  based  on  dis¬ 
tance  from  the  end  of  the  conuaction  ranging  from  2.91x106  to  5.%xl06.  Two  wall  inserts  were  used  to 
mount  the  microphones.  One  insert  held  both  a  single  microphone  and  the  surface-fence  wall-shear  stress 
gage  (which  will  be  described  later)  flush  with  the  wall  at  the  same  downstream  location.  It  is  shown 
with  the  surface  fence  in  place  in  Figure  3.  The  second  insert  was  the  same  insert  used  by  Blake  [1]  to 
make  convective  velocity  measurements.  It  is  shown  in  the  built-up  test  section  in  Figure  4.  It  consisted 
of  a  polished  circular  plug  flush-mounted  in  a  square  plate.  The  plug  had  a  line  of  1/8"  holes  along  its 
diameter.  By  placing  the  a  pair  of  microphones  in  selected  pairs  of  holes,  convective  speeds  and  coher¬ 
ence  of  pressure  waves  in  the  turbulence  could  be  measured.  During  all  these  measurements,  all  unused 
holes  were  filled  with  putty  and  made  flush  with  the  surface  by  shaving  off  with  a  razor  any  putty  which 
protruded  above  the  wall  surface.  The  microphones  were  brought  flush  with  the  insert  surface  by  sliding 
them  through  the  back  side  of  the  insert  until  they  protruded  above  the  tunnel  wall  surface  and  then  were 
pushed  back  to  the  wall  using  a  short  straight  edge  to  stop  the  motion  when  the  cap  surface  met  the  wall 
surface.  All  wall  pressure  measurements  were  made  with  Bruel  &  Kjaer  Model  4138  1/8"  microphones 
covered  by  pinhole  caps  having  a  hole  diameter  of  1/32". 

The  microphone  signals  were  amplified  by  a  pair  of  Bruel  &  Kjaer  Model  2607  measuring 
amplifiers  before  being  sent  out  of  the  wind  Uinnel  blockhouse  the  spectral  analyzer.  Before  each  run 
both  microphone/follower/amplifier/analyzer  channels  were  calibrated  using  a  Bruel  &  Kjaer  Model  4220 
pistonphone  which  produced  a  5(X)  Hz  tone  at  124  db.  Before  most  runs  the  microphones  were  baked 
for  forty-five  minutes  at  200oF  to  drive  off  any  moisture  which  might  be  present  on  the  diaphragm  sur¬ 
faces  and  cause  arcing. 

2.4.  Velocity  Measurements 

Profiles  of  the  streamwise  mean  and  fluctuating  velocity  fields  with  and  without  a  manipulating 
device  were  recorded.  A  TSI  Model  1261R  U-type  probe  placed  in  the  flow  was  balanced  by  a  DAN- 
TEC  Model  56C01/17  Constant  Temperature  Anemometer.  The  probe  was  mounted  on  a  DISA  Model 
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55A42  sting  and  inserted  through  the  side  of  the  tunnel  opposite  the  wall  pressure/wall  shear  measuring 
insert.  The  sting  position  was  precisely  set  by  mounting  it  in  a  screw  traverse  having  a  positioning  reso¬ 
lution  of  .005  inch.  It  is  shown  in  Figure  5.  The  anemometer  was  calibrated  against  a  pitot-static  probe 
placed  in  the  core  flow  of  the  tunnel.  The  pitot-static  probe  transducer,  a  Validyne  Model  CD23,  had  in 
turn  been  calibrated  against  a  Betz  water  manometer.  Both  the  Validyne  transducer  and  the  anemometer 
were  calibrated  before  each  profile  measurement.  The  anemometer  response  was  linearized  by  first 
recording  a  raw  volts  versus  dynamic  head  schedule  with  a  data  acquisition  program  on  the  HP- 1000 
computer  and  then  fitting  a  fifth  order  polynomial  to  the  curve.  All  test  measurements  were  translated  to 
velocity  using  this  curve. 

During  the  experiments,  velocity  sampling  was  done  with  the  HP- 1000  computer.  Approximately 
100,000  samples  were  taken  at  each  y  value  in  order  to  reduce  the  deviation  of  the  average  velocity 
from  le  true  velocity  to  .01  meters  per  second. 

2J.  Shear  Measurements 

The  measurement  of  mean  wall  shear  offers  some  keen  experimental  challenges.  The  primary 
challenges  have  to  do  with  intrusion  into  the  flow  field  and  calibration. 

A  sublayer  surface  fence  was  chosen  as  the  probe  for  wall  shear  measurement.  It  consists  of  a 
flush-mounting  plug  inside  of  which  is  a  plate  which  separates  two  chambers  and  which  protrudes  a  very 
short  distance  into  the  flow.  Each  chamber  has  a  port  to  a  pressure  transducer  and  another  port  to  the 
wall  surface  beside  the  fence.  When  a  boundary  layer  is  established  moving  perpendicular  to  the  fence, 
the  flow  separates  on  the  back  of  the  fence  and  a  pressure  drop  Ap  develops  across  the  fence.  If  the 
fence  is  short  enough  that  it  is  in  a  region  where  viscous  momentum  terms  dominate  the  convective 
terms,  the  response  of  the  probe  will  be  linear.  The  surface  fence  built  for  these  experiments  was 
approximately  the  height  of  the  viscous  sublayer  and  therefore  interacted  with  the  flow  field  in  a  way 
which  made  the  convective  terms  significant  The  calibration  curve  of  the  fence.  Figure  14  ,  shows  the 
wall  shear  stress  Tw<x:Ap2/3  ,  a  weak  non-linearity  for  all  but  the  smallest  values  of  Xw  .A  sketch  is 
shown  in  Figure  6.  It  has  been  described  by  Winter,  Femholz,  and  Head  and  Rechenberg  [14,15,16,17] 
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and  has  several  additional  advantages  over  other  measurement  techniques;  those  being  high  sensitivity 
and  directionality. 

The  benefit  of  remaining  inside  the  viscous  sublayer  is  that  the  measurement  no  longer  depends  on 
a  knowledge  of  the  shape  of  the  overlap  layer  but  rather,  as  Coles  [18]  has  pointed  out,  on  a  constant 
stress  layer  which  must  exist  very  near  the  wall  for  any  type  of  boundary  layer.  This  is  unlike  the  oft- 
used  Preston  tube  which  is  commonly  calibrated  in  flows  in  which  the  top  surface  of  the  tube  is  as  high 
as  several  hundred  y+  [19].  Such  a  calibration  then  requires  a  knowledge  of  the  behavior  of  the  "law  of 
the  wall"  and  how  it  might  be  modified  by  boundary  layer  manipulation. 

2.6.  Spectral  Mer  urements 

Point  spectra  were  measured  with  the  support  of  the  HP- 1000  computer  and  the  program 
P_SPECT4  which  contains  the  Fast  Fourier  Transform  (FFT)  subroutine  given  by  Newland  [20].  The 
subroutine  calculates 

One  hundred  spectra  each  having  2048  points  were  taken  and  averaged  at  each  flow  condition.  Data 
were  sampled  at  SO  kHz  and  the  resulting  bandwidth  was  12.2  Hz. 

Cross  spectra  were  measured  using  the  GenRad  Model  2515  16  channel  analyzer.  All  had  a  2  Hz 
bandwidth  and  1000  spectra  made  up  each  average.  The  continuous  cross-spectrum  is  expressed  as 


mdx 

The  analyzer  does  its  analysis  by  performing  a  direct  Fourier  transform  on  each  channel  and  then 

14.  Winter,  K.  G.,  (19T7),  "An  Outline  of  the  Techniques  Avail^le  for  the  Measurement  of  Skin  Friction  in  Turbulent 
Boundary  Layers".  Prog.  Aerospace  Sci.,  vol.  18.  pp.  1-57. 

15.  Dengel,  P.,  Femholz,  H.  H.,  and  Hess,  M.,  (1986),  "Skin-Friaioo  Measuremenu  in  Two  and  Three-Dimensionally 
Highly  Turbulent  Flows  WiUi  Separation",  Fint  Eutxtpean  Turbulence  Conference,  Lyon,  France,  (July  1986). 

16.  Vagt,  J.  D.,  and  Femholz,  H.,  (1973),  "Use  of  Surface  Fences  to  Measure  Wall  Shear  Stress  in  Three-dimensional 
Boundary  Layers",  Aeroruutical  Quarterly,  vol.  24.  part  2,  pp.  87-9t,  (May  1973). 

17.  Head,  M.  R.,  and  Rechenberg,  I.,  (1962),  "The  Preston  Tube  as  a  means  of  Measuring  Skin  Friction",  J.  Fluid 
Mechanics,  vol.  14,  pp.  1-17. 

1 8.  Coles,  D.  E.,  "The  Law  Of  The  Wall  In  Turbulent  Shear  Flow",  50  Jahre  Grenzschlchtforschung  ,  (H.  Gortler  &  W. 

Tollmien,  eds.),  F.  Vieweg  und  Sohn,  Braunschweig,  pp.  153-163. 

19.  Goldstein,  R.  J.,  ed..  Fluid  Mechanics  Measurements ,  Hemisphere,  New  York,  (1983),  pp.  559-611. 
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computing  the  auto-power  spectra  which  is  the  DFT  multiplied  by  its  complex  conjugate.  It  also  retains 
the  single  channel  DFTs  and  takes  complex  conjugate  products  to  get  the  phase  and  power  of  the  cross¬ 
spectra.  The  total  bandwidth  for  each  cross  spectra  was  approximately  twice  the  frequency  at  which  the 
coherence  of  the  signal  fell  below  .0025.  Both  signals  were  sent  through  anti-aliasing  filters  set  to  low 
pass  frequencies  below  the  sampling  rate.  All  spectra  were  taken  at  29  m/s. 


3.  EXPERIMENTAL  RESULTS 

3.1.  Mean  Velocity 

The  first  set  of  velocity  measurements  wctc  of  the  clean  (i.e.  unmanipulated)  tunnel  wall  boundary 
layer.  These  profiles  were  validated  in  much  the  same  way  Blake  [1]  did  in  his  smooth  wall  measure¬ 
ments.  The  dimensional  profile  was  plotted  in  velocity  versus  logy  coordinates  and  a  line  was  fitted  by 
regression  analysis  to  the  points  which  appeared  to  represent  the  "law  of  the  wall"  layer.  It  has  been 
mentioned  before  that  the  distance  of  the  probe  from  the  wall  could  be  measured  with  a  sufficiently  great 
precision.  It  was  not,  however,  sufficiently  "zeroed"  until  this  point  in  the  reduction.  Any  misjudgement 
in  the  zero  value  of  y  would  show  up  as  a  curvature  of  the  log  layer  on  the  semi-log  plot.  The  profile 
was  properly  "zeroed"  by  iteratively  adjusting  the  y  origin  until  the  most  colinear  set  of  points  resulted 
on  the  plot.  It  was  found  that  the  zero  value  on  the  traverse  changed  for  different  tunnel  speeds  and  the 
change  was  of  such  magnitude  and  direction  as  to  make  it  appear  the  wooden  walls  of  the  duct  were 
bending  slightly  and  progressively  inward  with  increasing  tunnel  speed. 

Figure  7  displays  the  boundary  layer  mean  profiles  nondimensionalized  by  outer-layer  variables 
and  is  a  demonstration  of  the  similarity  of  the  clean  wall  profiles  for  different  tunnel  Reynolds  numbers. 
Figure  8  shows  the  same  profiles  in  "law  of  the  wall"  coordinates  plotted  against  Coles  [18]  tabulated 
values  for  the  "standard"  mean  profile.  The  logarithmic  layer  is  described  by  the  equation 

20.  Newland,  D.  E.,  An  Introduction  to  Random  VlbraUons  and  Spectral  Analysis  ,  Longman,  Inc.,  New  York,  (1983), 
pg.  220. 
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The  value  of  friction  velocity  for  each  profile  was  determined  by  reguessing  the  "log  layer"  points  until 
the  value  ol  C2  resulting  from  the  regression  made  the  line  through  them  coincide  with  the  line  of  Coles. 

Figure  9  is  a  set  of  profiles  of  the  manipulated  boundary  layer  again  normalized  by  outer-layer 
variables.  In  all  cases,  the  shape  factor  has  increased;  typically  by  approximately  0.08  with  no  ticar 
trend  with  downstream  distance  being  apparent.  These  plou  are  now  not  necessarily  to  be  called  self¬ 
similar  because  new  length  scales  have  been  introduced  to  the  outer  variables;  those  being  primarily  the 
height  of  the  honeycomb  and  the  comb  diameter.  Figure  10  shows  the  same  profiles  scaled  by  comb 
height  What  is  apparent  from  the  plots  are  the  global  changes  in  tiie  profile  shapes  compared  to  the 
unmanipulated  ones.  What  is  apparent  from  the  data,  particularly  at  station  1,  is  the  "stepping"  of  the 
profile  left  by  the  internal  "boundary  layer"  growth  in  the  individual  combs.  The  profiles  at  station  2 
have  an  inflection  point  at  approximately  one  comb  height.  Near  the  wall,  the  value  of  ^  is  compara¬ 
tively  low  while  farther  from  the  wall,  near  the  outer  portion  of  the  boundary  layer,  it  is  comparatively 
high.  If  one  thinks  of  this  profile  as  an  indicator  of  mean  vorticity  in  the  boundary  layer,  the  manipu¬ 
lated  profile  has  shifted  the  vorticity  distribution  outward  from  the  wall  and  thereby  lessened  the  tran¬ 
sport  of  high  momentum  fluid  to  the  inner  portions  of  the  layer. 

Figure  11  is  a  set  of  "law  of  the  wall"  plots  of  the  manipulated  boundary  layer  at  various  down¬ 
stream  locations  plotted  along  with  the  Coles  data.  All  look  very  much  like  a  equilibrium  profile  which 
has  been  misscaled  by  using  a  friction  velocity  value  which  is  too  low.  All  of  the  data  are  shifted  in  the 
direction  of  larger  u+  for  the  same  y+  with  only  the  profiles  at  the  location  furthest  downstream  showing 
an  indication  of  strong  similarity  based  on  inner  variables.  The  profiles  taken  at  station  2  have  similar 
shape  but  their  separation  from  each  other  on  the  plot  in  such  a  regular  fashion  would  suggest  they  may 
scale  with  Reynolds  numbers  based  on  comb-related  variables  such  as  ma,.i(.aL.tcd  dispuicc.nent  thick¬ 
ness,  distance  downstream,  or  comb  height. 


The  skin-friction  coefficients  corresponding  to  these  friction  velocity  values  are  plotted  against 
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Reg.  and  Ax/5*  in  Figures  12  and  13  along  with  an  enipirical  relation  given  by  White  [21],  It  is 

Cf  =0.0128  Reg. -i 

White  estimates  the  curve  to  be  accurate  to  within  four  percent  and  the  data  are  at  most  three  percent 
from  the  line.  The  skin  friction  coefficients  of  the  manipulated  flow  are  from  24%  to  as  much  as  38%  of 
the  clean  flow  values  and  appear  to  have  the  trend  of  increasing  with  increasing  Reg..  Friction 
coefficients  appear,  for  the  region  measured  at  least,  to  be  decreasing  with  increasing  downstream  dis¬ 
tance.  It  would  be  expected,  however,  that  an  equilibnum  boundary  layer  would  reestablish  itself  by 
some  distance  downstream  and  that  the  friction  coefficient  would  return  to  the  expected  values. 

I  he  calibration  curve  for  the  surface  fence  is  shown  in  Figure  14.  It  shows  the  expected  and 
desired  linear  relationship  between  wall  shear  and  signal.  Again,  the  values  of  wall  shear  were  gotten  by 
fitting  of  the  respective  logarithmic  portions  of  the  profiles. 

3.2.  Turbulence  Levels 

Figures  15  and  16  are  plots  of  turbulertce  level  normalized  by  friction  velocity  versus  distance 
from  the  wall  normalized  by  displacement  thickness.  Figure  17  is  of  the  same  manipulated  profiles  but 
nondimensionalized  by  comb  height.  The  manipulated  profiles  have  comparatively  higher  intensities 
farther  from  the  wall  and  lower  levels  near  the  wall.  An  example  of  the  "plate"  effects  described  earlier 
is  that  just  downstream  of  the  comb  and  within  one  comb  height  of  the  wall,  the  manipulated  profiles 
have  a  certain  amount  of  "stagger"  to  them  because  of  the  individual  wakes  left  by  the  elements  of  the 
comb.  The  profiles  taken  at  station  2  have  the  expected  increase  coinciding  with  the  inflection  point  in 
the  profiles. 

3-3.  Wall  Pressure 

3J.1.  Point  Spectra 

Figure  18  is  a  plot  of  total  r.m.s.  pressure  at  certain  points  downtream  of  the  honeycomb  at  two 
tunnel  speeds.  The  dotted  and  dashed  lines  represent  r.m.s.  pressures  found  in  the  equilibrium  flow.  The 


21.  White,  F.  W.,  Viscous  Fluid  Flow  ,  McGraw-Hill,  New  York.  (1974),  pg.  495. 
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data  were  gotten  from  the  values  of  the  respective  autocorrelation  functions  at  zero  time  delay.  Because 
the  4138  Microphone  is  not  intended  to  measure  frequencies  below  10  Hz  and  because,  as  Blake  [1]  has 
reported,  large  amplitude  tunnel  disturbances  exist  below  70  Hz,  the  signals  were  high-passed  at  70  Hz 
before  being  analyzed.  The  maximum  reduction  in  wall  pressure  is  not  immediate  but  rather  approxi¬ 
mately  fifty  displacement  thicknesses  downstream.  In  air,  the  reduction  corresponds  to  as  much  as  4  db. 
Perhaps  there  is  a  zone  immediately  following  the  comb  in  which  the  wakes  shed  by  the  elements  of  the 
comb  dominate  the  turbulence  followed  by  a  zone  in  which  the  free  stream  velocity  and  pressure  gra¬ 
dient  begin  to  reestablish  an  equilibrium  flow. 

Pressure  spectral  densities  for  the  highest  three  tunnels  speeds  are  shown  in  Figures  19-22.  The 
spectra  were  smoothed  by  averaging  each  point  in  the  spectra  with  its  nearest  ten  neighboring  points. 
Every  tenth  point  is  plotted.  The  points  for  which  Corcos  [22J  has  predicted  a  3  db  reduction  in  trans¬ 
ducer  response  are  marked  with  a  caret  on  each  figure.  Nondimensionalization  is  done  in  the  first  three 
plots  with  free  stream  velocity,  mean  wall  shear,  and  displacement  thickness  and  in  the  second  with 
transducer  diameter  instead  of  displacement  thickness.  The  use  of  transducer  diameter  has  been  sug¬ 
gested  by  Foxwell  [23]  because  much  of  the  data  he  has  collected  for  a)5*/f/~  >1.3  scale  better  on  it 
than  on  displacement  thickness.  With  the  exception  of  station  1,  the  wall-shear  &  displacement  thickness 
nondimensionalization  makes  the  spectra  fall  very  nearly  one  atop  the  other.  Nondimensionalization  by 
transducer  diameter,  if  anything,  is  less  csq)able  of  scaling  high  frequency  data. 


3  J.2.  Streamwise  Coherence 

Figures  23-27  are  plots  of  streamwise  coherence  against  the  phase  of  the  cross  spectrum  along 
with  an  exponential  envelope  suggested  by  Farabee  for  equilibrium  flow.  Here  coherence  is  defined  as 


r(ri/3.(0) 


It  serves  as  a  measure  of  the  flow’s  ability  to  maintain  its  structure  as  it  convects  downstream.  When 


22,  Coroos,  G.  M.,  (1963),  "Resolution  of  P-essure  in  Turfjulence”,  J.  Acoustical  Society  of  America,  vol.  35,  no.  2,  pp. 
192-199. 

23.  Foxwell,  J.  H.,  (1966),  "The  Wall  Pressure  Specuum  Under  a  Turbulent  Boundary  Layer",  A.U.W.E.  Tech.  Note 
218/66,  (Aug.  1966). 
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plotted  against  phase  for  a  given  separation,  each  point  nominally  represents  the  phase  change  of  an 
eddy  of  a  given  frequency  as  it  convects  from  one  microphone  to  the  second  and  the  degree  to  which  it 
has  lost  Its  structure.  Phase  is  given  by 

a(o)o^  i)  =  tan-i 

WrftOoTl) 

In  all  cases,  there  is  a  more  rapid  decay  than  the  unmanipulated  flow  with  coherence  tending  to 
move  toward  equilibrium  with  increasing  downstream  distance.  It  is  possible,  once  again,  that  the  meas¬ 
urements  taken  at  station  1  are  in  the  zone  of  structures  which  are  of  the  order  of  the  comb  diameter  and 
which  typically  have  less  distance  to  travel  before  decaying. 

3  J  J.  Spanwise  Coherence 

Figures  28-30  are  similar  plots  of  spanwise  coherence  against  phase.  Here  no  convective  velocity 
can  be  discemable  from  the  measurements  and  Ucp  from  the  streamwise  cross-spectrum  for  the  same 
microphone  separation  at  the  same  manipulato'  position  is  used  along  with  (o  and  rs  to  calculate  a  phase 
relevant  to  the  convecting  structure.  The  notable  feature  of  these  plots  is  that  the  presence  of  the  mani¬ 
pulator  increases  spanwise  coherence  and  it  is  not  clear  why  this  is  true.  The  comb  may  be  shedding 
spanwise  vortices  off  its  top  over  a  broad  frequency  range. 

3.4.  Convective  Velocities 

Figures  31-35  show  pressure  wave  convective  speed  as  a  function  of  Strouhal  number  along  with 
the  line  reported  by  Bull  [24].  Phase  convective  speed  was  gotten  from  the  phase  component  of  the 
cross-spectrum  and  the  relation 

rr  _  OWi 

o^mTTT 

The  data  shows  little  difference  from  the  equilibrium  flow  until  (o8*/l/-  reaches  0.3,  the  number  at 
which  it  has  been  supposed  strucou'es  become  small  enough  to  convect  in  the  boundary  layer.  The  plot 
of  Ucp  for  the  largest  microphone  separation  is  treated  as  an  anomolous  case  and  may  be  the  result  of 


24.  Bull,  M.  K.,  Wilby,  J.  F.,  and  BUclonan,  D.  R.,  (1963),  "Wall  Pressure  Fluctuations  in  Boundary  Layer  Flow  and 
Response  of  Simple  Structures  to  Random  Pressure  Fields",  A.A.S.U.  Report  Vo.  2^3,  (July  1963). 
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the  microphones  being  placed  in  two  zones  of  the  flow  which  have  fundamental  differences.  That  the 
presence  of  the  manipulator  leaves  the  low  frequency  components  unchanged  may  tell  something  about 
their  source  location.  The  manipulator  dimensions  are  on  the  order  of  the  boundary  layer  thickness;  it  is 
an  outer  layer  device  and  it  would  be  expected  to  alter  most  or  all  of  the  outer  layer  structure.  The  low 
frequency  sources,  then,  may  be  in  the  inner  layers  and  associated  with  events  there,  such  as  bursting. 


4.  CONCLUSIONS 

1.  The  mean  wall  pressure  spectral  level  of  the  manipulated  flow  scales  better  with  mean  wall  shear 
stress  than  does  the  spectra  of  the  unmanipulated  flow. 

2.  The  manipulator  decreases  the  streamwise  and  increases  the  spanwise  coherence  above  (i>6*/f/»  = 
0.3. 

3.  The  manipulator  moves  the  active  region  farther  from  the  wall. 

4.  The  manipulator  reduces  convective  phase  velocities  for  frequencies  above  (o6*/t/-  =  0.3. 

5.  The  manipulator  reduces  mean  wall  shear  stress  and  root-mean-squared  pressure  for  downstream  dis¬ 
tances  of  the  order  of  100  displacement  thicknesses. 

6.  There  is  a  zone  immediately  downstream  of  the  manipulator  in  which  the  dominant  manipulator  length 
scale  is  the  comb  diameter  and  a  second  zone  further  downstream  which  is  dominated  by  the  comb 
height 
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Front  View 

Figure  la.  Honeycomb  Manipulator  Sketch. 
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Figure  3.  Point  Spectra  &  Shear  Measuring  Insert 


Figure  4.  Cross-Spectral  Insert  and  Test  Section. 
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Figure  5.  Hot  Wire  Traverse. 


Figxire  6.  Surface  Fence  Wall  Shear  Stress  Gage. 


Figure  10b.  u/V^  vs.  y//ic.  Manipulated  Flow.  Station  2. 
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40.  0  r  Figure  1  lb.  m+  vs.  y+  Manipulated  Flow.  Station  2. 
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Figure  15.  u'/u*  vs.  y/6*.  Equilibrium  Flow. 
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Figure  16b.  u'lu*  vs.  y/8*.  Manipulated  Flow.  Station  2. 
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Figure  16c.  u'fu*  vs.  y/6* .  Manipulated  Flow.  Station  3. 
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Figure  17b.  m7m*  vs.  y/Zic.  Manipulated  Flow.  Station  2. 


Figure  17c.  u'lu*  vs.  y//ic-  Manipulated  Flow.  Station  3. 
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Figure  20a.  Point  1 
^  Variables.  Station  1. 
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Figure  20b.  Point  1 
Variables.  Station  2. 
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Figure  22c.  Point  Pressure  Spectra.  Manipulated  Flow.  Transducer 
30. 0  Variables.  Station  3.  F 
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Figure  23.  Streamwise  Coherence  vs.  Phase,  ri/5*  =  .91 , 1.05 


Figure  24.  Streamwise  Coherence  vs.  Phase,  r  i/5*  =  2.19 , 2.53 


Figure  25.  Stream  wise  Coherence  vs.  Phase,  r  i/5*  =  4.01 , 4.63 


Figure  26.  Streamwise  Coherence  vs.  Phase,  r  i/5*  =  6.93 , 8.00 
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Spanwise  Coherence  vs.  Phase.  rs/S 


Figure  30.  Spanwise  Coherence  vs.  Phase,  r^/b*  =4.01 , 4.63 
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Figure  31.  Convective  Phase  Speed  vs.  Strouhal  Number,  ri/6 
1.05 


Figure  32.  Convective  Phase  Speed  vs.  Strouhal  Number.  r\/6*  =2.19 


Convective  Phase  Speed  vs.  Strouhal  Number,  r  i/6 


